Adiabaticity and localization in one-dimensional incommensurate lattices 
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We experimentally investigate the role of localization on the adiabaticity of loading a Bose-Einstein 
condensate into a one-dimensional optical potential comprised of a shallow primary lattice plus one 
or two perturbing lattice (s) of incommensurate period. We find that even a very weak perturbation 
causes dramatic changes in the momentum distribution and makes adiabatic loading of the combined 
lattice much more difficult than for a single period lattice. We interpret our results using a band 
structure model and the one-dimensional Gross-Pitaevskii equation. 



Disorder plays an important role in many condensed 
matter systems, [U |2 with deep connections to quantum 
chaos [3] , but can be difficult to systematically study due 
to the challenge of creating reproducible and quantifi- 
able disorder. The control available in ultra-cold atom 
systems [J makes it an attractive platform to study dis- 
order [5] m [3 IB . To date much of the work adding 
disorder to ultra-cold atom systems has explored time- 
independent properties, but the long timescales associ- 
ated with cold atoms allows investigation of dynamical 
properties as well (see [5[ |9[ UHl E] and ref. therein). 
In this work we examine the ability of a quasi-disordered 
system to adiabatically follow changes in the Hamilto- 
nian. The presence of disorder produces a complicated 
eigenvalue spectrum, which greatly affects the adiabatic- 
ity criteria. The physics of localization phenomena also 
has a significant impact on time-dependent processes, 
such as adiabaticity. Small perturbations to the Hamil- 
tonian can cause large changes to the ground state wave- 
function over large length scales, making it difficult for 
the system to adiabatically follow changes. One recent 
theoretical study shows that adiabaticity in gapless sys- 
tems is non-trivial, particularly in lower dimensions [12]. 
Here we show that even in a gapped system such as ours, 
adiabaticity is complicated by the presence of disorder. 

We study adiabaticity in a quasi-disordered system 
by adding one or two weak incommensurate lattices to 
a one-dimensional optical lattice loaded with a Bose- 
Einstein condensate (BEC). Localization occurs in both 
disordered and strictly incommensurate potentials [SI [131 
[Ml [T5] although with distinct differences, which tend to 
disappear in finite-sized systems such as ours. We ob- 
serve a complex momentum distribution of the atoms 
due to the presence of weak perturbing lattices follow- 
ing a ramped loading process that would be nearly adia- 
batic for a single lattice. We gain insight into the distri- 
butions from single-particle band structure, and observe 
that the effects of the perturbations disappear as inter- 
actions increase, as they suppress the long wavelength 
density modulation of the wavefunction. 

We form a BEC of -10"^ ^'^Rb atoms in the state 
\F = 2,m/ = 2) in a magnetic trap with ujx ^ 
^ 27rx410 Hz and co,. ^ 27rxl20 Hz. To reduce the 



effects of interactions, the trap is subsequently weakened 
giving final frequencies cOx ~ 27rx40 Hz, cj^y ^ 27rx20 Hz, 
and ujz ^ 27rx30 Hz. We load the BEC into a ID opti- 
cal potential, created by the addition of a primary and 
perturbing lattice(s). The total potential is 



Vtot = ^{^W + ^Iv' + ^l^') + Vi sin'(/^i^) 
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where M is the atomic mass, ki = 27r/A, and A = 796.6 
nm. The ratios k2/ki = 0.806 ± 0.002 and ks/ki = 
0.919 ±0.004 are extracted from images of atomic diffrac- 
tion. For the bulk of the experiments described here, the 
lattice depths are Vi = 4.6 ± 0.3 Erec and ^2,3/^1 = 
0.059 ± 0.003, respectively, where Free = n^kl/2M. 
To prevent interference between the lattices, we detune 
the beams several MHz from each other with acousto- 
optic modulators, so that the coupling terms between the 
beams oscillate rapidly compared to the atomic motion. 
The lattices are created by reflection off an in-vacuum 
gold mirror situated approximately 2 mm from the mag- 
netic trap center. Each beam is reflected at a different 
angle (l9i = 180°, 6>2 = 143°, and 0^^ = 155°) to yield 
a ID pseudo-disordered potential [16]. Because the gold 
mirror surface that defines the standing waves nodes is 
the same for all three lattices, they are phase- locked to- 
gether. We load the atoms into the combined lattices by 
ramping up the intensities linearly, keeping a fixed ratio 
between ¥2^3 and Vi. A ramp time of 1 ms is chosen for 
most experiments and is sufficiently long to ensure load- 
ing a magnetically trapped BEC into the lowest band of 
the primary lattice. In order to remain adiabatic, the ex- 
cited state population, |an,g(t)p, must stay <C 1, where 
n is the band index, and q is the quasimomentum. One 
can calculate a corresponding time scale using 
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which, for loading a ^^Rb BEC into a single lattice at 
ki ^ 27r/A, is satisfied for times much greater than 5 
/is [iTl HH mj. This short time scale is because the 
nearest excited band with allowed transitions is approxi- 
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mately 4 Ti^k^ /2M separated from the ground band with 
initial <7 = 0, even for an arbitrarily weak lattice with 
reciprocal lattice vector of 2k. As the lattice depth in- 
creases, the bands continue to separate. Adiabaticity be- 
comes difficult for q near the band edge, where the ini- 
tial energy gap vanishes. For our experiments, the initial 
momentum range is ±0.15 ki. Additionally, the single 
particle band picture begins to break down as the BEC 
becomes strongly interacting. 

To assess our ability to adiabatically load a lattice, we 
perform two types of experiments. In all experiments, we 
absorption image the cloud after 22 ms of time-of- flight. 
Following the ramped loading and a variable hold time, 
we abruptly turn off both the lattice and the magnetic 
trapping potentials. The image yields the momentum 
distribution of the atoms trapped in the lattice (interac- 
tions during time-of- flight are negligible). In the second 
method, we slowly ramp down the lattice and then release 
from the magnetic trap. If the process is fully adiabatic, 
the cloud should return to the momentum distribution of 
the the original BEC. For comparison, we present absorp- 
tion images after the sudden turn-off, with and without 
the perturbative lattice(s) in Fig. 1. The existence of 
the weak perturbing lattice (Fig.l(c-e)) markedly modi- 
fies the momentum distribution with the appearance of 
new peaks, even though the perturbing lattice is so weak 
that it does not produce any observable diffraction alone 
(Fig. 1(b)). In Ref. [13 it was suggested that an indica- 
tion of localization is the appearance of additional mo- 
mentum peaks in the matter-wave interference pattern. 

In order to understand why such a small perturbation 
has dramatic effects on the momentum structure and ex- 
cited fraction we use single-particle band structure. Even 
if the perturbing lattice is strictly incommensurate with 
the primary lattice {k2/ki = a, an irrational number), we 
can approximate a as a ratio of two large integers, f /g. 
In a finite system such as ours where the BEC occupies 
~ 70 sites of the primary lattice and does not extend 
over many periods of the beat frequency lattice (created 
by the combination of the lattice potentials) we expect 
this approach to yield reasonable predictions [20] . 

Assuming no excitations to higher bands, we predict 
the momentum distribution for each lattice configuration 
using a band structure calculation with f /g = 9/11 and 
f^g = 10/11, approximations to the experimental ratios. 
We present the results of this calculation over-layed with 
population amplitudes (green bars) extracted from fits of 
our data in Fig. 1. In the case of three lattices, experi- 
mental resolution limits our ability to resolve the struc- 
ture, and is represented by wider bars. The combination 
of the lattices creates a complex momentum structure, 
and the band structure calculation is in good qualitative 
agreement with our measurements. However, since the 
band-structure calculation assumes no excitations and a 
1 ms ramp causes depletion of the ground band, the pre- 



(I) 





Vl^4.6 Erec 






<> 


am- 



(11) 

0.6 
0.4 
0.2 




-2 2 

momentum (hk) 

FIG. 1: (I) Absorption images of a BEC loaded into different 
lattices for /j^bec ~ 100 Hz after 20 ms TOF: {ai)Vi = 4.6 Erec 
and (b) V2 = 0.3 Erec]{c) Vi + V2] (d) + + V3, with V2 
= V3 = 0.15 Erec ; (e) 7.5 Erec and 0.5 Erec lattice. We turn 
on lattices using a 1 ms linear ramp followed by a 1 ms hold 
and abrupt turn-off. (II) Using the data from the absorption 
images ((I)-(a,c-e)) we extract the relative population ampli- 
tudes and overlay them with a band structure calculation (no 
interactions) of the ground state momentum distribution. 



dieted populations do not quantitatively agree with the 
data for multiple lattices. This mismatch increases as 
the depths of both the primary and perturbing lattice(s) 
are increased. For deep lattices (Fig. l(i)), depending on 
the choice of //^ , band structure predicts peaks spaced 
at ki/g^ closer than the beat frequency, yet the envelope 
does not change. We find these disappear much more 
rapidly than the beat frequency peaks as interactions in- 
crease (see below). 

We also examine the changing wavefunction and the 
corresponding band-structure in Fig. 2. Using the 
/c— vector of the primary lattice to determine the size 
of the first Brillioun zone, we plot the first two bands 
for the case of one and two lattices. The dominant ef- 
fect of the perturbing lattice is to open small gaps in 
the ground band, as well as to slightly flatten the band. 
Starting at g = and traveling along the band, the 
first gap corresponds to quasimomentum of {ki — k2)- 
If instead one chooses the first Brillioun zone to span 
±(/ci — /C2), there exists a new band that is much lower 
in energy than the first excited band in the single lat- 
tice case, which determines a new energy scale relevant 
to adiabaticity. This new energy scale is also responsible 
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for localizing the wavefunction at the spatial period cor- 
responding to the beat frequency (Fig. 2(c) compared to 
Fig. 2(b)). In order to satisfy the adiabaticity condition 
(Eq. 2) for two lattices with our parameters {V2 ~ 0.059 
Vi, /ci//c2 ^ 0.8), the loading ramp must be much longer 
than 4 ms, a thousand times longer than for a single lat- 
tice. This is shown in Fig. 2(e), where we calculate the 
depletion of population in the ground band for different 
ramp rates as a function of primary lattice depth for the 
case of three lattices. The inset of Fig. 2(e) depicts an 
expanded view of ramp times yielding less than 5 percent 
excitation. Figure 2(f) shows the depletion of the ground 
state as a function of lattice depth for a ramp time of ap- 
proximately 2.5 ms for the case of one, two and three 
lattices. 

Since it is often the case experimentally that we start 
with a BEC with non-zero q (due to residual motion 
in the magnetic trap), we also show curves for loading 
a BEC with q = 0.1 ki into a potential comprised of 
two and three lattices. This value of q was chosen to 
lie within the new, smaller first Brillioun zone. These 
curves show that the effect of additional perturbing lat- 
tices is the dominant effect on adiabaticity criteria. For 
q = 0.1 ki^ the adiabaticity time scale approaches 10 ms, 
compared to virtually no change for a single lattice with 
initial q of 0.1 ki. We experimentally observe that by 
ramping up the potential and then ramping down with 
additional lattices, longer ramp times (5 ms compared 
to 250 /is) are qualitatively better (less excitations), but 
due to interactions, we are never able to be fully adi- 
abatic. Although excitation to higher bands may also 
produce new momentum peaks, for the data presented 
in Fig. 1(c), we predict that between 1 and 5 percent 
of the population is excited depending on initial quasi- 
momentum. Thus most of the observed changes in the 
momentum distribution can be ascribed to modification 
of the ground band wavefunction, i.e. localization. 

The effects seen here are a result of the incommensu- 
rability of the lattices. For experiments using commen- 
surate lattices [19 where two lattices had a carefully 
chosen ratio of 3:1, the criteria for adiabaticity with re- 
spect to the band structure is not significantly altered. 
In this case the higher order momentum peaks due to 
the perturbing lattice overlap the beat frequency peaks, 
which results in the wavefunction having modulation only 
at that spatial frequency. This difference can be seen by 
calculating the ground state wavefunction for the two 
cases. For approximately incommensurate lattices, the 
perturbation has a strong effect, with small spatial fre- 
quencies appearing (Fig. 2(c)). In order to further quan- 
tify this dramatic change, we plot the population of the 
p = peak, po, as a function of lattice depth for a sin- 
gle lattice, and for a primary lattice with a perturbative 
lattice (with V2 = 0.059 Vi). The black points are experi- 
mental measurements of the peak amplitude, which show 
a dramatic reduction of po when the perturbing lattice 
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FIG. 2: (a) Calculated first two bands for a single lattice and 
for two incommensurate lattices. Note the gaps which form 
in the ground band, (b-c) Plot of an example wavefunction 
multiplied by a Gaussian envelope (b) without and (c) with 
a perturbing lattice and (d) the probability amplitude, po, 
as a function of lattice depth for p = peak for the case of 
a single lattice and a combined potential. The changes in 
slope correspond to the depths at which the atoms become 
localized at different spatial frequencies, (e) Calculation of 
ground band population for differing ramp time ((i) 2.5 /is, 
(ii) 250 /is, (iii) 2.5 ms) for case of three lattices; (iv) for 
loading a single lattice with a ramp time of 2.5 ms. Inset in 
(e) is a zoomed-in view for comparison of (iii) and (iv). (f) 
Calculated ground state population for loading a (i) single 
lattice, (ii) two lattices, (iii) two lattices starting with ^ = 0.1 
and three lattices starting with (iv) ^ = and (v) ^ = 0.1 for 
a fixed ramp time of 2.5 ms. 



is added. We show two band structure calculations for 
//^=4/5 (dashed) and 9/11 (solid), both falling within 
the error bars of our experimental ratios. The signifi- 
cant difference between these cases is due to population 
growth in sub-beat frequency spaced momentum peaks 
predicted for the ratio 9/11. The error bars represent 
statistical uncertainties in populations and lattice depths 
(due to fluctuations in lattice beam intensity). We expect 
that because the turn-on of two lattices is not fully adi- 
abatic there will be additional depletion of the central 
peak due to excitations into higher bands, as observed. 

The production of low spatial frequency, long- 
wavelength components in the wavefunction with the ap- 
plication of weak incommensurate lattices yields insight 
into the connection between disorder, localization, and 
incommensurate systems. In the canonical view of local- 
ization, destructive interference due to reflections over 
long distances leads to the localization of the wavefunc- 
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tion [21 [21] . Here we see that a single incommensurate 
perturber yields a complex long wavelength structure. If 
we add a second perturber, getting closer to true disorder, 
we see even more momentum components. Experimen- 
tally, due to finite resolution, this appears as a broaden- 
ing of the distribution (Fig. 1(d)). Comparing three lat- 
tices to two, we observe a larger spread in the central fea- 
ture, as well as a corresponding decrease in optical depth. 
This indicates that there are unresolved peaks beneath 
the overall envelope. If something competes with these 
long range interferences, such as atom-atom interactions, 
we can expect localization phenomena to be drastically 
modified [11 . Indeed, until recently |5| [6] Anderson lo- 
calization had not been seen in cold atom systems due to 
the effect of atom- atom interactions. 

To study the effects of interactions in our system, 
we perform experiments in three different magnetic 
traps of differing frequencies with chemical potentials of 
/iB£;c/^^(2500,400, 100) Hz. We observe (Fig. 3) that 
as the overall confinement is increased, both the structure 
in the interference pattern such as the ki — k2 momentum 
peaks disappears and the overall size of the cloud, after 
TOF, increases. 

We compare the results of the non-interacting band- 
structure model to simulations of our system using the 
Gross-Pitaevskii equation (GPE). We reduce the 3D 
equation to ID and solve using the split-operator method 
in combination with imaginary time evolution [22 . The 
simulations using the GPE predict that in the moder- 
ately tight trap {jiBEc/h = 400 Hz), the ki — k2 (beat 
frequency) peaks slightly persist, however they disappear 
in the tightest trap. Fig. 3 shows absorption images 
with results of GP simulation in three different traps. In 
our system, we cannot decouple the size of the sample 
from interactions, which both contribute to the amount 
of spatial localization. This can be done if one utilizes 
a Feshbach resonance, as in [6 , to vary the strength of 
the scattering length. We calculate the effect of arbitrar- 
ily increasing the interactions in our weakest trap (where 
the structure is most apparent) and see that the sub-beat 
frequency peaks vanish rapidly, followed by the ki — k2 
peaks. This can be understood in that interactions drive 
transitions between states of different quasimomentum, 
washing out the discrete low momentum features, and 
destroying the long range spatial periodicity of the wave- 
function. 

In conclusion, we have presented data indicating that 
small perturbations to a one-dimensional lattice system 
in the form of quasi-disorder, while leading to localiza- 
tion of the wavefunction, also drastically changes the dy- 
namics of the system. We present a theoretical treat- 
ment which suggests that although we are only slightly 
modifying the energy of the system, the large alterations 
of the wavefunction demand time scales for adiabaticity 
that are orders of magnitude longer than for loading a 
single lattice. Because of the sensitivity of the adiabatic- 




FIG. 3: (I) (a-c) Absorption images in a traps with iibec ^ 
(100,400,2500) Hz with a primary plus one perturbing lat- 
tices. (II) (d-f) Results of GP simulation for (a-c) showing 
the effects of interactions. 

ity criteria in the presence of perturbations, one should 
be careful when studying disordered lattice systems to 
identify and characterize any forms of disorder present in 
the potential, intentional or otherwise. This work shows 
that disorder can have a strong influence on dynamics, 
and that the long timescales of cold-atom optical lattice 
systems makes them ideal for further explorations. 
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